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OPTIMUM  DETECTION  OF  RANDOM  SIGNAL  IN  NON-GAUSSIAN  NOISE 
FOR  LOW  INPUT  SIGNAL-TO-NOISE  RATIO 


INTRODUCTION 


A  finite-time  observation  is  made  of  a  weak  signal  (if  present)  in  a  non-Gaussian  noise 
background.  Optimum  detection  consists  of  computing  the  likelihood  ratio  (LR)  of  the  observed 
waveform  samples  and  comparing  the  LR  with  a  fixed  threshold,  which  enables  a  statement 
about  signal  presence  versus  absence  in  that  time  interval.  The  random  noise  samples  are 
denoted  by  {n*}  for  k  =  1  :K.  It  is  presumed  that  these  noise  samples  are  statistically 

independent  with  common  first-order  probability  density  function  (PDF)  pn  (x).  If  the  signal  is 

present,  its  samples  are  {sk } . 

For  a  deterministic  signal  {sk },  the  LR  and  log  LR  take  the  forms 


LR  =  Y\  ~  ^  logLR  =  ]T  [tf(x*  -J*)-«(xt)], 

k=]  Pn\Xk )  k-l 


(1) 


where  {x* }  is  the  received  data,  and  q(x)  =  log  p„  (x).  For  low  input  signal-to-noise  ratio  (SNR), 
the  approximation 


logLRs-£  »'(*»)  ** 


k= 1 


f  KMS 
h  p„m 


(2) 


follows.  Thus,  for  low  SNR,  the  dominant  term  in  the  optimum  processor  depends  on  the 
derivative  of  the  q  function  and  is  linear  with  the  signal  samples  {sk}. 


On  the  other  hand,  for  a  stationary  random  signal  with  zero  mean ,  the  optimum  processor  for 
low  input  SNR  must  depend  on  higher-order  statistics  of  the  input  signal.  At  the  very  least,  the 
covariance  of  the  input  signal  will  be  required  or  must  be  estimated.  Also,  it  can  be  expected 
that  higher-order  information  about  the  non-Gaussian  noise  PDF  pn  (x)  will  be  required  in 
addition  to  q'(x).  This  is  indeed  the  case.  To  make  reliable  decisions  about  signal  presence  or 
absence,  these  more  stringent  information  requirements  for  the  random  signal  case  require 
additional  data  on  both  the  received  signal  and  the  noise. 
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DERIVATION  OF  LIKELIHOOD  RATIO  PROCESSOR 


Samples  {x^ }  for  k  -  1  :K  are  available.  Under  the  hypotheses 
H0:xk  =nk, 

Hl:xk=  sk+nk,  ^ 

the  optimum  detector  of  signal  presence  versus  absence  is  desired.  Noise  sequence  {n* }  consists 
of  independent  samples,  with  known  first-order  (FO)  PDF  pn  (x),  which  is  arbitrary.  Random 
signal  sequence  {sk }  is  characterized  by  joint  ^-dimensional  PDF  ps  (w, , . . . ,  uK ).  The  optimum 
LR  detector  forms  the  ratio 


LR  =  )  _  Psn(Xl>‘’'>Xx) 

k=\ 

for  comparison  with  a  fixed  threshold. 

The  numerator  of  equation  (4)  is  given  by  A'-dimensional  convolution: 

psn(xl,...,xK)=  J-  J  dux--duK  ps(ux,...,uK)  pn{xx -Ux)-p„  (xK -uK).  (5) 

For  low  input  SNR  and  zero-mean  processes,  the  joint  PDF  ps  (m,  , . . . ,  uK )  is  concentrated  near 
the  origin  of  ^-dimensional  u-space,  relative  to  the  spread  of  the  remaining  FO  noise  PDFs 
iPn(xk  ~uk))-  Therefore,  a  typical  noise-PDF  term  in  equation  (5)  can  be  expanded  near  uk  =  0 
according  to 

PMk~uk)  =  pn(xk)~»k  tiM  +  Wk  P'Mk)  for  k  =  l:K.  (6) 

This  expansion  is  to  second  order  in  uk. 


Two  auxiliary  functions  are  defined  as 


PM 

Pn(*y 


b(x)  = 


PM 

Pn(x)' 


Then,  equation  (6)  becomes 


(7) 


2 


(8) 


pn(xk-uk)  =  pn (xk )[l  +  uk  a(xk)  +  \u2k  b(xk )]  for  k  =  \:K. 

Substitution  of  this  expansion  in  equation  (5)  and  retention  of  terms  through  second  order  in 
{uk },  which  is  consistent  with  second-order  expansion  (6),  yields 


P^.^TliPnM}  f  "  | dux--duK  ps(u{,...,uK ) 


i+2>*  «(**)+tZm*  6(x*)+iZw*  uj  a(xk)a(*j) 

k  k  k*j 


(9) 


The  three  types  of  integrals  in  equation  (9)  can  be  evaluated  as 

J—  \dul—duK  ps(ul,...,uK)uk  =  jduk  psl(uk)uk  =0, 

J-- J dux—duK  ps(ux,...,uK)u2k  =  J duk  psX(uk)u2k  =  a2,  (10) 

J--- J dux--duK  ps(uv...,uK)uk  Uj  =  Jj \duk  duj  ps2(uk,Uj)uk  Uj=a]  pkJ, 

where  a2  is  the  common  signal  variance  of  samples  {sk},  and  pkj  is  the  normalized  signal 
covariance  coefficient  of  sk  and  Sj. 

The  use  of  equation  (10)  in  equation  (9),  followed  by  substitution  in  equation  (4),  leads  to 
the  approximation 

LRsl  +  Jof  2*(xt)+ioJ  2  Pa  a(xk)a(*j)  forlow  inPut  SNR-  (n) 

k  k*j 

Since  data-independent  additive  constants  and  positive  scalars  do  not  affect  the  receiver 
operating  characteristics  (ROCs)  of  a  processor,  an  equivalent  processor  to  equation  (1 1)  is 

K  K 

. ^)  =  ^(h)+I^  a(xk)a(Xj),  (12) 

*=1  k,j= 1 

where  modified  covariance  matrix  p  =  [pkj  ]  is  obtained  from  K  x  K  matrix  p  -  [pkJ  ]  by  replacing 

all  the  unity  diagonal  elements  with  zeros.  This  modification  enables  equation  (12)  to  be 
expressed  in  the  compact  form 

K 

L(xx,...,xK)  =  ^b(xk)  +  aT  pa,  (13) 

k=\ 


3 


where  K  x  1  vector  a  -  [a(x1)--a(xAr)]r.  This  is  the  (approximate)  optimum  detector  for  low 
input  SNR  and  any  joint  signal  statistics.  Although  the  signal  level  is  not  involved  in  equation 
(13),  the  normalized  signal  covariance  matrix  p  must  be  known  to  employ  equation  (13). 

An  alternative  form  for  the  optimum  processor  is  obtained  by  defining 

q(x)  =  logp„(x).  (14) 

Then,  using  equation  (7), 


»’<*) .  PM  .  _aW>  q.(x)  _  PMriM-  W  m  b(x)_a(xy 

Pn  \X)  P„(X) 

The  optimum  processor  in  equation  (12)  can  now  be  expressed  as 

K  K  K 

L(x„...,xK)  =  YjbM-Yja(xkf  +  £  pkj  a(xk)a(Xj) 

*=1  *=!  k,j=\ 

£ 

=  Y,q'M  +  E  Pkj  q'Mq'iXj). 


(15) 


(16) 


An  arbitrary  constant  can  be  added  to  q\x )  without  affecting  the  performance  of  this 

As  a  special  case,  if  the  signal  samples  {sj  are  uncorrelated,  then  matrix  p  =  0, 
equation  (13)  reduces  to 


4(x 


1»* 


•^K)  =  ^b{xk)  = 


*=1 


y1  7>/"(x^) 

hpn(*kj 


processor. 

and 


(17) 


Nonlinearities  q’{x)  and  a(x)  are  not  involved  at  all  in  this  processor.  However,  more 
generally,  when  the  signal  samples  are  correlated,  then  matrix  p  &  0  and  the  last  term  in 
equations  (12)  and  (13)  allows  for  interaction  between  the  {x* }  terms  by  means  of  the  nonlinear 
transformation  a(x)  -  -  p'n(x)lpn(x). 


The  noise  samples  {nt }  have  been  taken  statistically  independent  in  the  derivation  above. 

If  the  sampling  rate  used  to  generate  the  data  is  increased,  the  noise  samples  will  become 
statistically  dependent  on  each  other.  This  effect  is  investigated  in  appendix  A  for  a  simple 
example.  It  is  shown  that  the  deflection  criterion  is  insignificantly  improved  through  an 
increased  sampling  rate.  Sampling  at  the  highest  rate  that  still  results  in  uncorrelated  samples 
will  extract  all  the  detection  capability  that  a  given  data  segment  length  contains.  Higher 
sampling  rates  just  require  more  data  processing,  with  little  if  any  improvement  in  performance. 
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GAUSSIAN  NOISE  PDF 


pM) 


1  xJ 

=—  exp  -  — 


P'n  0)  =  ~^Pn  (X),  Pn  (X)  =  (  ~~T - T  I Pn  (*)• 

< 7 -  °n  °n  J 


a(x )  =  ,  b{x )  =  \  (x2  -  (j\  ), 


4  Zfc  -°"«)+X  a,  xkxt  . 

L  k  k*J 


Discarding  additive  constants  and  positive  scalars,  an  equivalent  processor  is 

2X + e  Pkj  xk  x7  =  E  pkj  x*  xj  =  *T p  x> 


where  X  =  [x,  •  •  •  xK  ]T .  This  is  a  “generalized”  energy  detector,  weighted  by  the  normalized 
signal  covariance  matrix/?. 

Alternatively,  from  equations  (14)  through  (16), 

q(x)  =  -irT -log(V2?r<T„),  q'(x)  =  -^Y,  q\x)  =  -\, 

2(7 »  °»  (20) 

K  i  K 

i(Xj , . .  )  —  2~  ^  Pkj  5 

O'* 

which  is  tantamount  to  processor  (19).  The  gff(x)  term  has  effect  on  the  performance  of  the 
processor  for  the  case  of  Gaussian  noise. 

The  results  in  equations  (19)  and  (20)  hold  for  low  input  SNR,  g]Ig\,  and  any  signal  joint 
PDF  ps  («j  ,...,uK  ).  It  is  informative  to  compare  these  results  with  the  exact  LR  processor  for 
Gaussian  signal  and  Gaussian  noise.  The  latter  processor  is,  from  equation  (4)  and  discarding  a 
positive  scalar, 


LR  =  exp  --Xr*-,X  +  -!TXrX  =exp  ±-XT\\l-Rj;  X 
2  2<t„  2  l  <r„  J 


where  signal-plus-noise  covariance  matrix 


(22) 


For  low  input  power  SNR^/cr*  « 1,  the  inverse  of  matrix  Rsn  is 


Substituting  in  equation  (21)  yields 


LR  =  exp 


This  is  simply  a  monotonic  transformation  of  equivalent  processors  (19)  and  (20);  this 
equivalence  holds  under  approximation  (23). 

More  precisely,  using  the  identity 

I  -  (I  +  rp)~ 1  =rp(I  +  rp)~' ,  r  =  a)  I c r2„  <  1, 

equation  (21)  becomes  exactly 

XTp(I  +  rpy'X  . 


An  equivalent  exact  LR  processor  is,  for  this  Gaussian  noise  PDF, 

L(xl,...,xK)=XTp(I  +  rpyiX  =  XTp  ( I-rp  +  r2p 2 - )X 

=  XT p  X-rXr p2  X  +  r2  XT p3  X - . 


(23) 


(24) 


(25) 


(26) 


(27) 


Terms  beyond  the  first  involve  additional  factors  containing  powers  of  input  SNR  r  =  <r2Jcrl 

and  matrix/?,  both  of  which  allow  for  dropping  these  terms  for  low  input  SNR  with  negligible 
effect  on  detection  performance. 

Expansion  (9)  was  conducted  only  through  second  order.  In  appendix  B,  the  third-  and 
fourth-order  expansions  are  derived  for  those  interested  in  pursuing  these  higher-order  terms 
numerically. 
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EXPONENTIAL  NOISE  PDF 


/>„(*)  =  7^exp  -^1  °l=1P\ 

^P  V  P ) 

?W  =  -^-l°g(24 

P  (28) 

qXx)  =  {Z  &rrI<00}  =  -/gn<X)’  fmaJ & 

q\x)  =  -jS(x). 

The  difficulty  in  this  example  is  due  to  noise  PDF pn(x)  not  possessing  a  derivative  atx  =  0.  The 
last  term  in  equation  (16)  becomes  — r^.Pu  sgnfx^.)  sgn(xy ),  which  does  not  appear  to  have  an 

P  kj 

immediate  problem.  However,  the  first  term  involving  q\x)  obviously  needs  special  treatment 
and  is  taken  up  in  appendix  C. 


CAUCHY  NOISE  PDF 

B  1  x  P  -2jc  .  p  6x2  -ip1 

Pn(*)  = - T—tf’  Pn(x)  =  -T-2 >  Pn(x)  =  -J~ y 

ft  x  +P  ft[x2+fi2)  x{x2+fi2) 


2x  l/„\  _  6*2  -  'ip2  ^_2(*1 ~P4 


a(x)=  j  2-,  b(x)  =  — - f-,  q\x)  =  -a(x),  q\x) 


x2+P 1 


(S+/32)1 


(*’+/>T 


Large  xt  values  are  suppressed  by  both  components,  and  q",  of  the  optimum  processor: 

?'(**)  — — .  q\xk)~\- 


ALTERNATIVE  NOISE  PDF  (< a  =  p'3  42 1  ft  ) 

pUx)=a4^j’  p'Ax)~a4x2v^r 


4x3  t/  N  „25x4-3 p*  „  ,  ,  2  x4-3/?4 

"(x)=7^’ tw=4x  «<x)=4x 
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Again,  large  xk  values  are  suppressed  in  a  similar  fashion  to  equation  (30).  But,  small  xk  values 
are  also  suppressed  by  a(x),  which  has  cubic  behavior  near  the  origin.  Odd  function  a{x)  has  a 
peak  at  x  =  3U4/3  of  value33/4/^. 

Some  additional  impulsive  noise  models  are  given  in  appendix  D. 


MULTIPATH  EXAMPLE:  GENERAL  NOISE  STATISTICS 

Letting  the  received  signal  contain  a  direct  path  component  and  a  delayed  (surface  bounce) 
multipath  component: 

sk  =d*  ^>0,  (32) 

where  zero-mean  sequence  {d*}  consists  of  uncorrelated  samples, 

E{dkdj}  =  (r2d  S(k,j).  (33) 


Then,  signal  covariance 

E{*kSj}  =  (J]  [8{k,j)  +  a  8(k-m,j)  +  a  8(j -m,k)  +  a2  8(k-m,j -«)], 
leading  to  total  received  signal  power  a 2  =  a]  (1  +  a2)  and 


1  0  •••  0  (3  0  ••• 

0  10  0/3  0 

i  0  \  \ 

P  =  |0  ,  /3  = 

P  0 

0  (3  \ 

:  0 


a 


1  +  a" 


(34) 


(35) 


where  value/?  =  a /(l  +  a2)  appears  in  the  super-  and  sub-diagonals  that  are  m  removed  from  the 
main  diagonal.  Optimum  processor  (12)  or  (13)  now  takes  the  form 


K  2/y 

Z(x, ,...,  x* )  =  YJb{\k )  +  -—f  £  a(xk )  a(xk+m ) 
^=i  1  +  a 


(36) 


for  low  input  SNR,  where  nonlinear  transformations  a(x)  and  b(x)  are  given  by  equation  (7). 
The  first  term  in  equation  (36)  is  a  sum  of  samples  {xk }  subject  to  nonlinearity  b(x),  while  the 
second  term  is  an  autocorrelation  (at  delay  m )  of  distorted  samples  {a(xk )}.  Realization  of 
processor  (36)  requires  knowledge  of  scaling  a  and  delay  m  in  model  (32).  Lack  of  this 
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information  will  require  several  guesses  of  both  parameters  and  repeated  evaluation  of  the 
second  term  in  equation  (36)  to  realize  a  maximum  ofI(x1,...,xA.).  However,  the  latter  sum  in 
equation  (36)  does  not,  by  itself,  require  knowledge  of  a,  while  its  scale  factor  depends  solely 
on  a.  This  observation  can  be  used  to  simplify  the  computation  of  equation  (36)  for  multiple 
guesses  of#  andm. 

An  alternative  expression  for  the  optimum  processor  in  equation  (36)  is  available  upon  use 
of  equation  (15);  namely, 


k  k  2a 

L(xx , . . . ,  xK )  =  %q\xk )  +  £?'(x*  )2  + 

k= i  k= i  i  + a 


K-m 

Y,q\xk)q\x 

k+m  )* 


k= 1 


If  signal  delay  information  m  and  signal  strength  information  a  are  unknown,  a  suboptimum 
processor  is  available  upon  dropping  the  last  term: 

K  K 

L(xlt...,xK  )  =  £y(x*)  +  ^q\xk)2. 

t= i 
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ESTIMATION  OF  OPTIMUM  NONLINEARITIES  FROM  DATA 


To  implement  the  optimum  processor  in  equations  (13)  or  (16),  it  is  necessary  to  determine 
the  nonlinearities  q'(x )  and  q"{x).  When  an  analytic  expression  is  available  for  noise 

PDF  pn  (x),  as  in  equations  (18),  (28),  (29),  and  (31),  the  appropriate  derivatives  can  be  taken  and 

used  in  equation  (15).  This  information  is  needed  in  equation  (16)  in  addition  to  the  signal 
covariance  matrix  p. 

However,  when  an  analytic  noise  PDF  is  not  available,  an  estimate  of  the  noise  PDF  must  be 
made  from  a  finite  length  of  noise-alone  measured  data  (n^ }.  Furthermore,  two  derivatives  of 

the  noise  PDF  must  be  extracted  somehow  from  this  finite  data  without  exacerbating  the  effect  of 
taking  differences  of  noisy  estimates.  The  method  adopted  here  is  to  obtain  a  histogram  from  the 
measured  data,  fit  it  by  means  of  a  spline  with  several  continuous  derivatives,  and  then  take  two 
derivatives  of  the  spline  fit.  In  particular,  the  “spaps”  procedure  of  MATLAB  was  used  for  all 
the  following  examples. 

As  a  first  example,  100,000  independent  normalized  Gaussian  RVs  were  generated  and 
subjected  to  the  spaps  procedure.  The  results  for  the  three  functions  q  =  q(x),  q{  =  q'(x),  and 

^2  =  q”(x)  are  shown  in  figure  1.  An  estimate  for  each  function  is  displayed,  as  well  as  the  exact 
nonlinearity  from  equation  (20).  The  agreement  among  the  three  pairs  is  excellent  over  the  full 
range  of  ±4.25  standard  deviations  plotted.  The  corresponding  results  for  the  smaller  sample  size 
of  10,000  data  points  is  given  in  figure  2;  the  deviations  between  the  exact  and  estimated  results 
are  now  obvious  but  not  excessive.  Still  smaller  sample  sizes  lead  to  considerably  worsened 
estimates  of  the  three  q  functions. 

The  second  example  of  interest  is  a  non-Gaussian  RV  obtained  by  means  of  the  nonlinear 
transformation 

y=x+ax\  (37) 

where  x  is  a  Gaussian  RV.  Scale  factor  a  is  set  equal  to  0.1 .  For  a  unit- variance  RV  x,  the  non- 
Gaussian  RV  y  has  the  exact  q  functions  displayed  in  figure  3;  the  derivations  as  well  as  a 
program  for  these  functions  are  presented  in  appendix  E.  The  heavier  tails  of  this  PDF,  relative 
to  the  Gaussian  case,  are  evident  in  the  upper  q  function  in  figure  3,  which  was  plotted  over  the 
range  of  ±5  standard  deviations.  The  q 2  function  has  a  pronounced  dip  of  substantial  value  at  the 
origin,  while  the  qx  function  behaves  like  a  limiter  for  input  amplitudes  greater  than  unity. 
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Figure  2.  q  Functions  for  10,000  Gaussian  Data  Points 
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Figures  4,  5,  and  6  illustrate  the  type  of  stability  and  accuracy  that  can  be  attained  from 
independent  sample  sizes  of  le7,  le6,  and  le5,  respectively,  when  the  PDF  and  q  functions  must 
be  estimated  from  finite  data  segments.  The  results  in  figure  4  for  1  e7  samples  agree  very  well 
with  the  exact  results  in  figure  3,  except  on  the  edges  near  ±5  standard  deviations.  This  edge 
effect  is  to  be  expected  because  of  the  rarity  of  encountering  these  values  and  the  concurrent 
instability  of  the  histogram  estimate. 

The  results  in  figures  5  and  6  display  much  poorer  agreement  between  estimation  and  truth, 
even  though  the  two  upper  q  functions  appear  to  be  fairly  well  fitted  by  the  spline  procedure. 

The  problem  lies  in  getting  the  higher-order  derivatives,  which  accentuate  the  minor  deviations 
in  the  fits.  This  non-Gaussian  example  illustrates  that  a  considerable  amount  of  data  is  going  to 
be  required  if  reliable  estimates  of  the  qx  function,  and  especially  of  the  q2  function,  are  to  be 
obtained  directly  from  the  data.  This  is  not  an  unexpected  result;  estimation  of  higher-order 
derivatives  from  noisy  data  will  always  be  difficult  to  accomplish  successfully  unless  a  very 
significant  amount  of  independent  data  is  available  and  is  carefully  processed. 

Comparison  of  figures  1  and  6  also  reveals  that  estimation  of  the  q  functions  for  the  non- 
Gaussian  noise  case  is  considerably  less  stable;  equivalent  accuracy  in  the  estimated  results  for 
the  non-Gaussian  noise  case  will  require  more  data  than  for  the  Gaussian  noise  case. 
Furthermore,  lack  of  sufficient  data  may  preclude  decent  estimation  of  the  q2  function,  with  the 
attendant  degradation  of  performance  of  the  resultant  suboptimum  processor.  Observe  that  for 
the  Gaussian  noise  case,  the  q2  function  is  a  constant  and  can  be  totally  ignored  in  the  optimum 
processor.  This  is  not  the  case  for  any  non-Gaussian  noise  situation,  where  the  q2  contribution 
can  be  significant  in  making  the  proper  detection  decision. 
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PERFORMANCE  OF  SEVERAL  NONLINEAR  PROCESSORS 


Suppose  that  RV  x,  which  has  FO  PDF  px(x),  is  nonlinearly  transformed  by  monotonically 

increasing  function/,  yielding  RVy  =  / (x).  Letting  /  denote  the  inverse  function  of/  the  FO 
cumulative  distribution  function  (CDF)  of  RV  y  is  given  by 

=  Pr(y  <y)  =  Pr(/(x)  <y)  =  Pr(x  <  f(y))  =  cx(f{y)).  (38) 

The  corresponding  FO  PDF  is  given  by 

Py(y)  =  ^;cy  (y)  =  ^cx(f{y))  =  px(f{y))  f\y).  (39) 

But,  since 

y  =  f{x)  =  f(f{y)),  1  =  f\f{y))f\y),  (40) 

the  end  result  for  the  FO  PDF  of  RV  y  follows  as 


py{y) 


pXhy) ) 
f'Giy)) ' 


(41) 


This  last  step  avoids  the  need  to  evaluate  the  derivative  of  inverse  function  /,  as  in  equation 
(39),  which  is  often  a  much  more  complicated  nonlinear  function  than / 

As  an  example,  consider  the  cubic  nonlinear  transformation 

y  =  f(x)  =  x  +  ax\  a>  0;  f\x)  =  1  +  3ax2.  (42) 

Then,  the  inverse  of  this  transformation  is 


The  positive  square  root  is  taken  to  ensure  increasing  monotonicity  with  y.  Once  PDF  px  (x)  is 
specified,  all  the  quantities  needed  in  equation  (41)  are  available. 
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The  detailed  solution  of  the  cubic  nonlinear  transformation  in  equation  (42)  proceeds  as 
follows:  let  x-t  +  u,  where  t  and  u  are  arbitrary.  Substitution  in  equation  (42)  yields 


(t  +  u)(\  +  3atu)  +  a(t3  +  u3)  =  y. 

Now  let  u  =  -1/(3 at),  in  which  case  equation  (45)  reduces  to 
1 


(45) 


at3  - 


-  y,  t6  -—t3 - ^—7-  =  0, 

27  a  21a 3 


u  = 


2a  y 
-1 


\2  Qj 


1 


( 


21a3 


t  = 


y 


2  a  y 


ry? 

V2  aj 


xl/3 


1 


21a 


(46) 


3  at 


,  x  —  t- 


3  at 


This  is  the  result  quoted  in  equations  (43)  and  (44). 

The  optimum  nonlinearities  corresponding  to  PDF  py  (y)  in  equation  (41)  are 

q(y)  =  log py(y),  ql(y)  =  q'(y),  q2{y)  =  q\y)-  (47> 

Analytic  expressions  for  these  quantities  are  complicated.  Appendix  E  provides  a  MATLAB 
program  using  symbolic  math  for  these  tasks.  Simpler  analytic  fits  to  qx(y)  and  q2 (y)  are  also 
included  for  purposes  of  significantly  reducing  the  numerical  evaluation  of  equation  (47).  Figure 
7  depicts  the  optimum  nonlinearities  and  their  fits  for  an  example  with  a  =  0.1  and  a  Gaussian 
RV  x  with  zero  mean  and  unity  standard  deviation. 

A  comparison  was  made  of  the  ROCs  for  several  processors,  using  K  =  1000  independent 
data  samples  and  various  combinations  of  nonlinearities  qx  and  q2  in  figures  8  through  11.  For 
these  simulation  results,  1  million  trials  were  used.  Processor  A  used  only  the  linear  term 
qx=v  with  q 2  =  0,  where  r  is  the  received  data;  processor  B  used  qx  and  q2  equal  to  the 
optimum  fits  above;  processor  C  used  qx  as  the  optimum  fit  and  q2  =  0;  and  processor  D  used 
qx  =0  and  q2  equal  to  the  optimum  fit.  To  achieve  false  alarm  probability  Pf=  le-3  and 
detection  probability  Pd  =  0.5,  the  required  input  SNRs  for  processors  A  through  D  were  -10.1, 

-13.3,  -13.1,  and  -12.4  dB,  respectively.  Instead,  to  achieve  operating  point  (le-4, 0.9),  the 
required  input  SNRs  for  the  four  processors  were  -7.8,  -10.8,  -10.5,  and  -9.9  dB,  respectively. 
Finally,  for  processors  B,  C,  and  D,  the  required  input  SNRs  at  operating  point  (le-5, 0.99)  were 
-9.2,  -8.9,  and  -8.4  dB,  respectively.  In  all  cases,  processor  B,  which  used  optimum  fits 
qx  and  q2 ,  performed  the  best,  although  processor  C  trailed  by  only  0.3  dB;  that  is,  ignoring  the 
contribution  of  nonlinearity  q2  was  not  very  detrimental.  However,  processor  D,  which  uses 
only  q2 ,  was  more  costly  in  terms  of  loss  of  SNR.  It  should  also  be  observed  that  processor  A, 
which  used  only  a  purely  linear  term,  suffered  a  significant  3-dB  loss  relative  to  the  best, 
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processor  B.  Accordingly,  it  is  very  important  to  include  the  “saturating  effect”  of  the  first  type 
of  optimum  nonlinearity  qx ,  which  tends  to  suppress  the  large  outliers  of  the  cubed  Gaussian 
process  in  equation  (42). 
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SUMMARY 


The  optimum  processor  for  a  weak  stationary  random  signal  in  independent  non-Gaussian 
noise  has  been  derived  to  second  order  in  the  signal  strength.  The  information  required  to  realize 
the  likelihood  ratio  processor  consists  of  the  first-order  probability  density  function  of  the  noise 
process  and  the  covariance  function  of  the  signal  process.  More  precisely,  the  noise  density  and 
its  first  two  derivatives  are  required. 

Estimation  of  these  two  required  derivatives  of  the  noise  density,  from  a  finite  segment  of 
noise-only  data,  is  not  a  simple  task.  The  approach  adopted  here  is  to  fit  the  noise  histogram  by 
means  of  a  spline  with  several  continuous  derivatives.  Then,  analytic  derivatives  are  taken  of 
this  fit.  For  a  Gaussian  noise  process,  the  second  derivative  need  not  be  utilized,  making  this  a 
simpler  special  case.  For  a  non-Gaussian  noise  process,  both  derivatives  of  the  noise  density 
must  be  estimated,  making  the  estimation  task  considerably  more  difficult. 

Attempts  were  made  at  dropping  the  second-derivative  term  in  the  optimum  processor 
because  of  its  difficulty  of  reliable  estimation.  Simulations  of  the  receiver  operating 
characteristics  for  the  corresponding  suboptimum  processors  revealed  a  loss  of  several  decibels 
in  detectability  performance.  Accordingly,  to  achieve  optimality  it  is  necessary  to  conduct  the 
lengthy  data  collection  to  enable  reliable  estimation  of  the  higher-order  noise  density  derivatives 
and  characteristics. 
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APPENDIX  A 

EFFECT  OF  SAMPLING  INCREMENT  ON  PERFORMANCE 


The  effect  of  the  sampling  increment  on  the  deflection  criterion  in  the  detection  of  a 
deterministic  signal  in  colored  noise  is  of  interest.  The  received  waveform,  for  signal  present,  is 
s(t)  +  n(t),  where  s(t)  is  the  known  deterministic  signal  waveform  and  n(t)  is  a  stationary  noise 
process.  Data  processing  consists  of  sampling  this  received  waveform  and  cross  correlating  it 
against  the  stored  reference  according  to 

k= 1 

where  K  is  the  total  number  of  samples  utilized.  For  zero-mean  noise,  the  mean  value  of  RV  y 
is,  for  signal  present, 


My)  =  Z^(4)2-  (A~2) 

k 

The  variance  of  y  is,  under  both  hypotheses, 

c-!U)  =  l4) *fj ) K (h  - 1, )>  (A-3> 

kj 


where  Rn  (r)  is  the  correlation  function  of  the  stationary  noise  process  n(t).  The  quantity  of 
interest  here  is  the  dimensionless  deflection  criterion 


d  = 


My) 

o-(y)’ 


(A-4) 


In  the  special  case  of  a  Gaussian  noise  process,  this  statistic  completely  governs  the  ROCs;  more 
generally,  it  is  a  useful,  simple  measure  of  detectability  performance. 

For  equispaced  samples  with  increment  A,  tk  =  kA,  equation  (A-4)  becomes 


X  s(kA)2 

_ k _ 

X s(kA)s(jA)Rn((k-j )A) 


(A-5) 


The  signal  waveform  will  be  taken  to  be  a  constant  s0  over  the  fixed  observation  interval  of 
duration  T.  Then,  sampling  increment  A  =  T/K  and  equation  (A-5)  simplifies  to 


A-l 


For  K  2  FT ,Q  —  1 .  This  corresponds  to  A  =  T/K  =  1  /(2 F),  which  is  just  the  sampling  increment 
that  results  in  uncorrelated  samples  of  this  white  noise  process  with  normalized  correlation 
sinc(2F t).  This  can  be  considered  to  be  the  standard  sampling  increment  for  this  particular 
process.  For  other  values  of  K,  the  summation  in  equation  (A-8)  must  be  carried  out  in 
numerical  detail;  the  results  are  presented  in  figure  A-l  for  2FT=  100.  It  will  be  observed  that 
increasing  K  above  2FT=  100  does  not  significantly  improve  the  quality  ratio.  In  fact,  K  can  be 
decreased  to  about  60  before  a  significant  loss  is  observed. 

Figure  A-2  displays  the  corresponding  quality  ratio  for  a  noise  process  with  a  Gaussian- 
shaped  spectrum.  Again,  sampling  faster  than  necessary  does  not  increase  the  detectability  of  the 
signal  in  noise;  the  quality  ratio  simply  saturates  very  quickly  once  the  sampling  increment 
realizes  approximately  uncorrelated  samples.  Any  higher  sampling  rate  simply  requires  more 
data  processing  with  no  attendant  improvement  in  detectability  performance. 

It  should  be  observed  that  no  assumptions  have  been  made  about  the  statistics  of  the  noise 
process  n(f),  such  as  a  Gaussian  PDF.  The  conclusions  above  apply  to  stationary  noise  with 
arbitrary  first-order  statistics.  It  is  expected  that  similar  behavior  holds  for  weak  random  signals 
(rather  than  deterministic  signals)  in  noise,  when  correlators  or  energy  detectors  are  used  for 
processing. 
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More  generally,  as  K  tends  to  infinity,  the  second  factor  in  equation  (A-6)  tends  to 


For  very  large  FT,  this  quantity  approaches 


TRM 
3.(0)  . 


(A-9) 


(A- 10) 


where  Gn  (/)  is  the  noise  spectrum.  For  the  Gaussian-shaped  spectrum  in  figure  A-2,  the 
normalized  correlation  function  is  pn  (r)  =  exp(-4;r  F2t2)  and  the  spectrum  is 


(  n  f2  \ 

Gn  (/)  =  G„  (0)  exp  -  —  .  At  frequencies  /  =  ±F,  the  relative  spectral  value  is  down  by 
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exp(-;z74)  =  0.456  =  -3.41  dB.  The  second  factor  in  equation  (A-6)  again  tends  to  (2 FT)112  for 
large  FT,  as  K  increases.  The  reason  that  the  curves  in  figures  A-l  and  A-2  go  just  above  value  1 
is  due  to  the  factor  (1  -  |r|  /  T)  in  equation  (A-9),  which  slightly  attenuates  the  integral  in  the 

denominator;  the  limit  of  1  is  approached  only  for  very  large  FT. 
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APPENDIX  B 

HIGHER-ORDER  TERMS  IN  EXPANSION  (9) 


To  determine  the  third-  and  fourth-order  terms  in  expansion  (9)  in  the  main  body  of  this 
report,  the  bracketed  expansion  (8)  is  extended  to 


1  +  Uk  a(?-k  )  +  2  b(*-k  )  +  6  'Ujc  C&k  )  ■*"  24  Uk  k  ^ 


(B-l) 


where  auxiliary  functions 

P„(x)  P„(x) 

The  product  of  the  K  terms  of  equation  (B-l)  contains  the  third-order  quantities 


(B-2) 


r Zc(x*) a(xy)  uk  UJ  »  IXx')  a^)  a(*j)  ui  uk  uj  . 

k  k*j  l<k<j 


(B-3) 


and  the  fourth-order  quantities 


k  k*j  k<j 

(B-4) 

k*l 

i2>(x,)  u]  J>(x*)  a(*j)  uk  Uj ,  Za(x<»)  a(xt)a(xk)  a(\j)  um  u ,  uk  Uj . 

I  k<  j  m<l<k<j 

When  the  quantities  in  equation  (B-3)  are  substituted  in  equation  (9)  in  the  main  report,  the 
averages  over  the  joint  signal  PDF  ps(ul,...,uK)  will  require  knowledge  of  the  general  third- 

order  signal  moment  E(s,  sk  s y )  for  all  l,k,j.  If  the  signal  samples  in  equation  (3)  in  the  main 
body  of  this  report  are  taken  from  a  zero-mean  Gaussian  process,  these  third-order  moments  will 
all  be  zero. 

On  the  other  hand,  the  averages  over  the  quantities  in  equation  (B-4)  will  require  knowledge 
ofthe  general  fourth-order  signal  moment  E(sm  s,  sk  sy)  for  all  m,l,k,j.  For  a  zero-mean 

Gaussian  signal  process,  this  fourth-order  moment  can  be  expressed  as  combinations  of  second- 
order  signal  covariances  coefficients,  namely, 

<*s  {Pml  Pkj  Pmk  Plj  +  P mj  Plk  )’  ^ 

For  other  processes,  the  fourth-order  moment  is  much  more  difficult  to  evaluate. 
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APPENDIX  C 

MODIFICATION  OF  APPROXIMATE  LIKELIHOOD  RATIO  PROCESSOR 
FOR  EXPONENTIAL  NOISE  DENSITY 


When  the  additive  noise  in  model  (3)  in  the  main  body  of  this  report  has  an  exponential 
PDF,  equation  (28)  indicates  that  the  q\x)  component  of  the  optimum  processor  (16)  becomes  a 
delta  function.  This  unacceptable  term  is  due  to  expansion  (6)  breaking  down  in  the 
neighborhood  of  x*  =  0,  where  the  exponential  PDF  does  not  have  a  derivative.  A  closer  look 
at  equation  (5)  reveals  that,  for  low  input  SNR,  the  desired  operation  is  a  local  smoothing  of  the 
joint  noise  PDF  in  the  neighborhood  of  the  point  x„...,xr  This  suggests  that  nonlinearities 
q'(x)  and  q\x)  in  equation  (28)  be  smoothed  in  x  by  a  function  of  width  approximately  that  of 
the  standard  deviation  <rs  of  the  signal.  Since  the  smoothing  operation  in  equation  (5)  would 
tend  to  retain  the  shape  of  the  noise  PDF,  the  following  candidate  replacement  functions  were 
employed,  respectively,  instead  of  q\x)  and  q"(x) : 

qx  (x)  =  -  [1  -  exp(-|  x  |  /g)]  sgn(x),  q2  (x)  =  -  -  exp(-  |  x  |/g) .  (C-l) 

& 

These  functions  are,  respectively,  odd  and  even,  with  q2  (x)  being  the  derivative  of  qx  (x), 
consistent  with  that  same  property  of  original  nonlinearities  q'{x )  and  q”(x).  The  positive  scalar 
g  is  taken  as  c  <rs ,  where  c  is  of  the  order  of  1 .  Numerical  simulations  of  optimum  processor 

(16),  with  replacements  (C-l),  yielded  a  significant  improvement  in  the  ROCs  compared  with 
those  for  the  q\x)  nonlinearity  in  equation  (28);  the  q"(x)  nonlinearity  in  equation  (28)  had  to 
be  completely  discarded  because  it  cannot  be  implemented.  A  gain  of  about  4  dB  was  obtained 
by  use  of  replacements  (C-l)  in  optimum  processor  (16)  using  only  q'(x).  An  additional 

experiment,  in  which  the  q2  (x)  term  was  discarded,  led  to  a  degraded  performance  in  the  ROCs; 
accordingly,  it  is  concluded  that  the  q2  (x)  term  is  a  significant  contributor  to  the  approximate 
optimum  processor  and  must  be  retained.  The  ROC  performance  using  replacements  (C-l)  was 
not  very  sensitive  to  the  exact  value  employed  for  scalar  c;  thus,  exact  knowledge  of  the  input 
signal  standard  deviation  <js  is  not  a  necessity.  Also,  alternative  nonlinearities  to  those  in 
equation  (C-l)  were  tried,  without  significant  changes  or  degradation,  provided  the  proper 
scaling  g  was  employed. 

This  replacement  procedure  is  an  ad  hoc  approach  based  upon  smoothing  the  original  direct 
results  in  equation  (28).  To  lend  credence  to  this  approach,  an  exact  derivation  of  the  LR 
processor  (4)  for  a  Gaussian  signal  and  exponential  noise  was  conducted  for  K  =  2.  The  exact 
statistics  employed  are 
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i  [Ml 

Pn{'x)  =  2p\py  0’”"2>92’ 


2)  = 


2n<j*(\-p2)1 


M|2  +«2  -2pu{u 


2  r  1/9 1**  2 

2/1  ^2\ 


2^0-/»  ) 


(C-2) 


The  exact  LR  is  given  by  equations  (4)  and  (5)  with  K  =  2.  To  evaluate  this  double  integral,  the 
following  two  integral  results  are  very  useful: 


x2  +y2  -2 pxy 


(C-3) 


00  co 

J<ftJ</yexp(- jax2  -\Py2  +yxy  +  px  +  vy) 


2n  fay2  +  ftp1  +2ypv\  T . 

=7dexp — i^HI(u’v’p)’ 


(C-4) 


where 


i  =  aP-r\u  =  ]i{a-e^\v  =  ]^b-!^\P  =  pL 


Definition  (C-3)  is  given  in  reference  1,  section  26.3.  The  result  in  equation  (C-4)  is  obtained  by 
making  linear  shifts  of  the  variables  x  and  y  to  eliminate  the  linear  variables  in  the  integrand  and 
then  scaling  to  manipulate  the  double  integral  into  the  form  of  equation  (C-3). 

Upon  substituting  equation  (C-2)  into  equation  (5)  from  the  main  body  of  this  report, 
splitting  the  double  integral  into  four  parts  corresponding  to  the  noise  joint-PDF  breakpoints  at 
its  origin,  and  using  equation  (C-4)  four  times,  the  following  result  for  the  exact  LR  is  obtained 
for  as  >0  (after  considerable  manipulations): 


r  =  y’y.  =^’y2  =^>Hy.Hy2k  =rQ  +  p),tm  =r(l-p), 
e,  =  exp(r  +  a  -  y ,  -y2),e2  =exp (rtm  +a  +  y,  -y2), 

e3  =  exp(/-  tp+a  +  y]  +y2),e4  =exp(r/„,  +a-y,  +y2),  (C-5) 

L,  =  L(tp  -y,  fr,tp  -y2 /r,p\  L2  =  L(tm  +y,  /r,tm  -y2  lr-p\ 

L3  =  L{tp  +y, /r,tp  +y2  /r,p),  L4  =  L(tm  -y,  /r,tm+  y2  !r-p\ 

T-^(X| ,x2)  L|  ~i"  e2  L2  -t*  e3  L3  e4  L4 . 
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A  computer  routine  for  the  joint  Gaussian  exceedance  probability  distribution  Lih^k^p)  enables 
rapid  numerical  evaluation  of  this  LR  for  any  xt ,  x2 . 

The  approximate  LR  processor,  obtained  from  equation  (16)  by  use  of  replacements  (C-l),  is 
given  by 

K  K 

LRq(xl,...,xK)s'£q2(xk)  +  YjPkj<l\  (x*)4i(x,)-  (C'6) 

A=1  k,j=  1 


Several  comparisons  between  exact  result  L/?(x, , x2 )  and  LRq  (x, , x2 )  are  made  in  figures  C-l 
through  C-4  for  various  values  of  parameters  r,pn,c,xx/<Jn,x2/crn.  Since  additive  constants 
and  scaling  factors  do  not  affect  the  ROCs  of  a  processor,  LRq  (x, ,  x2 )  was  transformed  so  that  it 
had  the  same  values  as  ZJ?(x, ,  x2 )  at  the  origin  and  +  oo.  The  close  similarity  of  the  two  sets  of 
numerical  results  indicates  that  the  modified  procedure,  obtained  by  smoothing  the  nonlinearities 
q\x)  and  q\x)  when  necessary,  is  an  appropriate  and  viable  procedure;  this  is  especially 
pertinent  when  K  is  large  and  the  AT-fold  integration  of  equation  (5)  is  intractable.  Then, 
modification  (C-6)  becomes  an  attractive  and  useful  alternative. 

It  was  noted  in  equation  (17)  in  the  main  body  of  this  report  that  the  optimum  processor 
conducts  a  single  sum  of  the  transformed  input  data  {xk }  when  the  input  signal  samples  {st }  are 

uncorrelated.  However,  for  exponential  input  noise,  equation  (28)  reveals  that  the  pertinent 
nonlinear  transformation  involves  a  negative  delta  function,  which  is  meaningless.  Instead,  an 
exact  analysis  for  a  Gaussian  signal  PDF  yields  LR 


K 

KxIf...,xr)=n 


k=\ 


exp(jr2  +|y*|  +  y*)<D(-r-y*/r) 
+  exp(|r 2  +  |y  *  |  -  y  * )  ®(-r  +  y  Jr) 


(C-7) 


where  O(w)  =  jdt  exp(-t2/2)/V2^  is  the  normalized  Gaussian  cumulative  distribution 

—oo 

function.  Plots  of  the  k-th  term  on  the  right-hand  side  of  equation  (C-7)  versus  xk  reveal  a 
function  that  is  a  continuous  approximation  to  the  negative  delta  function  at  the  origin,  with  a 
width  proportional  to  signal  standard  deviation  crs.  This  is  further  confirmation  of  the  smoothing 
procedure  that  is  recommended  earlier  in  this  appendix  for  nonlinearities  that  are  meaningless  or 
cannot  be  realized  in  practice.  This  procedure  can  be  useful  for  other  noise  PDFs  pn  (x)  with 
non-analytic  behavior  at  some  argument(s). 
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APPENDIX  D 

TWO  IMPULSIVE  NOISE  MODELS 


The  first  noise  model  is  the  square  of  an  exponential  RV,  followed  by  a  multiplication  by  an 
equally  likely  ±1 .  This  RV  can  be  generated  by  the  operation 


(\  V  C24)1/4 
u  =  ±  -log(r)  ,  b  =  — — , 
b  y  o 


(D-l) 


where  RV  r  is  uniformly  distributed  over  (0,1),  and  cr  is  the  desired  standard  deviation  of  RV  u. 
The  PDF  of  RV  u  is 


b  exp 

p(u )  = - 7= -  for  all  u, 

4  l\u\ 


(D-2) 


and  its  even  moments  are  given  by 


E{  u2"}  = 


(D-3) 


The  odd  moments  of  u  are  zero. 


The  EDF  of  u  is  given  by 


e(u)  = 


j  exp (-byfu)  for  u  >  0 
1  -  \  exp(-6^jw| )  for  u  <  0 


(D-4) 


The  characteristic  function  of  u  is 


/(£)  =  f  du  exp(i^M)  p(u)  =  V ~n  Rej  ^  r==~  w  I  for  £  *  0. 


(D-5) 


This  is  a  closed  form  in  terms  of  the  w(z)  function  (see  reference  1,  chapter  7). 
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The  second  noise  model  is  the  fourth  power  of  a  zero-mean  Gaussian  RV,  followed  by  a 
multiplication  by  an  equally  likely  ±1.  Letting  q  be  a  ±1  RV  with  equal  probabilities,  for  an 
arbitrary  RV  x,  the  RV  y  =  q  x  has  the  following  properties: 

/>(£  I  q)  =  £,  {exp(/£qx)}  =  fx  (4q), 

Py  (u)=i  Px  («)  +  t  Px  (-«)  for  all  w,  ^ 

ey(u)  =  2ex(u)  +  \cx(~u)  for  all  w. 


For  the  particular  case  of  RV  x  =  g4,  where  g  is  a  zero-mean  unit-variance  Gaussian  RV,  the 
following  functions  govern  the  statistics  of  RVs  x  and  y: 


/*(£)=  f  du-jL=exp(-±u2  +  ifr*)=-l=LL=exp 
J  4ln  32£ 


f  :  ^ 


v32^y 


(  ;  A 


K 


1/4 


V32^y 


for  £  >  0, 


/,«)  =  Re 


l«expUk 


(  ,  ^  /  ,  A 


K, 


1/4 


v32l4y 


for  £  *  0,  /JO)  =  1, 


1 


Px(u)  =  —^<t>(u  )  for  w  >  0;  ex{u)  =  2  0(-«l/4)  forw>0, 

(M)  =  !3/4  ^(H,/4)=  ,3/4  exP(“  2  JH)  for  all  u  *  0, 


ev(u)  = 


8|wf -  '  V2^8|W| 

<t>(-wl/4)  fori/>0 


3/4 


[<D(|w|1/4)  for  w  <  Oj 
The  variance  of  RV  y  is  105.  All  the  odd  moments  are  zero.  The  functions 
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<f>{x)  =  {ln)An  exp(-x2 12),  <D(x)  =  fdt  </>(t), 


(D-8) 


are,  respectively,  the  PDF  and  cumulative  distribution  function  of  a  normalized  Gaussian  RV. 


D-2 


APPENDIX  E 

PROGRAM  FOR  OPTIMUM  q  FUNCTIONS 


clear  all  %  Optimum  ql (y)  and  q2(y) 

syms  y  a  pi  %  for  Y  =  X  +  a  XA3 

sx=  1 ;  %  Standard  deviation  of  X 


tl=y/(2*a); 

b=l/(3*a); 

t=(tl  +sqrt(tl  A2+bA3))A(l/3); 

x=simplify(t-b/t); 

fp=l+3*a*xA2; 

px=exp(-.5  *  (x/sx)  A2)/sqrt(2  *pi*  sxA2) ; 

py=simplify(px/fp); 

q=log(py); 

ql=diff(q,y); 

ql=simplify(ql); 

qli=inline(ql); 

q2=diff(ql,y); 

q2=simplify(q2); 

q2i=inline(q2); 

ad=.  1; 

yd=[0:.02:10]'; 

qld=qli(ad,yd); 

figure(l),  elf 

plot(yd,qld,'r') 

grid  on,  hold  on 

q2d=q2i(ad,yd); 

plot(yd,q2d,'b') 

yd=[0:.2:10]'; 

fld=-l  .6*yd./(l+.5*yd.Al  .6); 
plot(yd,fld,'r.') 


%  Solution  of  y  =  f(x)  =  x  +  a  xA3 
%  f  (x) 

%  Gauss  PDF(x)  of  X 
%  PDF(y)  of  Y  =  f(X) 

%  Optimum  ql(y) 

%qli(a,y) 

%  Optimum  q2(y) 

%  q2i(a,y) 

%a  =  0.1  example 


%  Fit  to  ql(y)  for  ad=.l 
%  Fit  to  q2(y)  for  ad=.l 


f2d=1.6*((yd/1.79).A2-l)./(l+1.6*yd  A2+.2*yd  A4); 

plot(yd,f2d,'b.') 

xlabel('y') 
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